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Abstract 

We have investigated the glassy behavior of a molecular crystal built up with chloroadaman- 
tane molecules. For a simple model of this molecule and a rigid fee lattice a MD simulation 
was performed from which we obtained the dynamical orientational correlators S\y (q, t) 

and the "self" correlators Sy),(t), with A = (£,m), X' = {£' ,wl). Our investigations are 
for the diagonal correlators A = A'. Since the lattice constant decreases with decreasing 
temperature which leads to an increase of the steric hindrance of the molecules, we find a 
strong slowing down of the relaxation. It has a high sensitivity on A, A'. For most (£, m), 
there is a two-step relaxation process, but practically not for (£, m) = (2, 1), (3, 2), (4, 1) 
and (4, 3). Our results are consistent with the a-relaxation scaling laws predicted by mode 
coupling theory from which we deduce the glass transition temperature T C MZ) = 217 K. 
From a first principle solution of the mode coupling equations we find T C MCT = 2Q7K. 
Furthermore mode coupling theory reproduces the absence of a two-step relaxation process 
for {£,m) = (2, 1), (3,2), (4, 1) and (4,3), but underestimates the critical nonergodicity 
parameters by about 50 per cent for all other {£, to). It is suggested that this underestimation 
originates from the anisotropic crystal field which is not accounted for by mode coupling 
theory. Our results also imply that phonons have no essential influence on the long time 
relaxation. 
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1 Introduction 



The mode coupling theory (MCT) [1] is the most successful microscopic approach 
to describe glassy dynamics of supercooled liquids on a qualitative, and partly even 
on a quantitative level [2,3,4]. Particularly, first-principle comparisons between the 
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solutions of MCT equations and results from experiments and simulations for vari- 
ous systems like the binary Lennard- Jones liquid [5], liquids of diatomic molecules 
[6,7], water [8,9], silica melt [10] and a model for orthoterphenyl [11] have con- 
firmed the quality of MCT. These comparisons were restricted to the glass order 
parameters, the nonergodicity parameters. Going beyond that, even time dependent 
quantities, like the intermediate scattering function, have been compared with each 
other for the binary Lennard- Jones liquid [12], binary mixtures of hard spheres 
[13,14] and the polydisperse quasi-hard-sphere system [15], again demonstrating 
consistency with MCT. 

Glassy dynamics of systems with self-generated disorder is not restricted to liq- 
uids. Molecular crystals in their plastic phase exhibit glassy behavior as well. This 
has been found experimentally three decades ago [16]. Recently two of the present 
authors have extended MCT to molecular crystals [17]. Using the static structure 
factors from Percus-Yevick theory [18] the orientational glass transition of uniax- 
ial hard ellipsoids on a simple cubic lattice has been investigated [17]. There it 
has been found that this transition is not driven by an orientational cage effect, 
analogous to supercooled liquids, but by the growth of orientational order. This 
growth manifests itself in an increase of the static orientational correlators at the 
Brillouin center or edge, which in turn leads to an increase of the memory kernel. 
However, before the corresponding orientational correlation length diverges at the 
corresponding equilibrium phase transition line, the nonlinear feedback mechanism 
of MCT results in an orientational glass transition. The same mechanism has been 
identified for a liquid of uniaxial hard ellipsoids [19]. Accordingly, the MCT glass 
transition line is located within the orientationally disordered ergodic phase [17,19] 
and not in the supercooled regime. 

One may ask whether there also exists a cage-effect-driven glass transition for 
molecular crystals. Good candidates are plastic crystals which undergo a first or- 
der equilibrium phase transition to an orientationally ordered phase. Several such 
systems exist such as cyanoadamantane, ethanol, cyclooctanol, difluorotetrachloro- 
ethane or Ceo [20]. Another one is chloroadamantane, which exhibits an equilib- 
rium phase transition at T e exp = 2AAK [21]. A simple model for chloroadamantane 
(see next section) has been studied by MD simulations [22,23]. A MCT analysis 
for the critical amplitudes and the a-relaxation time have demonstrated consis- 
tency with MCT predictions [23]. These predictions can really be extended from 
liquids to plastic crystals because the MCT equations for the latter have the same 
mathematical structure than for multi-component simple liquids [17]. 

It is the main goal of the present paper to perform a first principle comparison 
between the results for chloroadamantane from a MD-Simulation and MCT. The 
outline of our contribution is as follows. In the second section we will describe the 
model, introduce the orientational correlation functions, give some details of the 
simulational procedure and shortly describe the relevant MCT equations. Results 
are presented and discussed in the third section and the final section contains a 
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summary and some conclusions. 



2 Model and technical details 



Chloroadamantane C\qH 15 CI is a rather huge molecule which belongs to the sub- 
stituted adamantane family. It shows a plastic phase structure isomorphous to cyano- 
adamantane, but the chloroadamantane molecule possesses a smaller substitute and 
a faster dynamics well adapted for MD simulation investigations. Chloroadaman- 
tane undergoes at T ~ 2AAK a first order transition from an ordered monoclinic 
structure to a rotator phase with face-centered-cubic (fee) symmetry [21]. The 
plastic-liquid transition occurs at T m ~ AA2K [21]. 



Table 1 

Parameters for the two-site chloroadamantane model. 



Site-Site 


P 


q 


e (kJ/mol) 


(7(A) 


CI -CI 


12 


6 


1.441 


3.350 


CI - Adm 


14 


8 


3.087 


4.786 


Adm - Adm 


16 


11 


12.47 


6.200 



The simulated system is composed of rigid linear molecules with two sites: one 
chlorine atom (noted CI) and one super atom (noted Adm) that models the adaman- 
tane part Ci Hi 5 . The moment of inertia is 302.733 amu.A 2 . Molecular dynamics 
calculations were performed on a system of iV = 256 molecules (4 x 4 x 4 fee 
crystalline cells) interacting through a Lennard- Jones short range site-site potential 
of the form 

v{r) = Ae((a/r) p - (a/r) q ) 

where r is the distance between two different sites. The parameters e, a, p and q are 
specified in table 1 . 

The chloroadamantane molecule possesses a relatively large dipolar moment fl 
(2.39 Debyes) which is parallel to the molecular axis. The electrostatic interac- 
tions were handled by the Ewald method with two partial charges (q = ±0.151e) 
localized on both sites. Newton's equations of motion were solved with a time 
step of At = 5 fs. We worked in the NPT (constant number of molecules, tem- 
perature and pressure) and NVT (constant number of molecules, temperature and 
volume) statistical ensembles with periodic boundaries conditions. The sample was 
first equilibrated in the NPT ensemble. Then, MD runs at constant volume using the 
average volume determined from the NPT simulations were performed. 
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Fig. 1. Temperature dependence of the lattice constant obtained from experiments (full 
circles) and MD-simulation (open circles) 

Since our molecule-model has one rotational symmetry axis we can use the angles 
Qn(t) = (Q n (t), (j> n (t)) to characterize the orientation of the n-th molecule at site 
n and time t. The microscopic local orientational density is given by: 

Pn (n,t) = 5(n\n n (t)) (i) 

with 5(Q\Q') = (sin 0) _1 5(0 — Q')5((j) — 0'). Expansion with respect to spherical 
harmonics Y x (9, 0), A = (I, m) and performing a lattice Fourier transform leads 
to the tensorial orientational density modes: 

N 

p A (q,t)=*< Yl Y x (Q n (t))e^ (2) 

n=l 

where iV is the number of lattice sites and R n the lattice vector of the nth lattice 
site. The reader should note that we assume a rigid lattice and that q is restricted to 
the 1. Brillouin zone. Introducing the fluctuations 5p\(q, t) = p\(q, t) — (p\(q_, t)) 
we can define the time dependent, tensorial orientational correlators: 

47T 

Sxxi^t) = —(5p* x (q,t)5py(q,0)) (3) 

where (...) denotes the canonical average over the initial conditions. Note that the 
absence of a head-tail symmetry leads to nontrivial correlators for all £ > 1, while 
the absence of phonons leads to S\\>(q, t) = for A = (0, 0) and/or A' = (0, 0) 
within the 1 . Brillouin zone. 

These correlators have been investigated by a MD-simulation. In contrast to ear- 
lier work [22,23] the present simulation has been done for each temperature T for 
a rigid lattice with lattice constant a(T) determined from the average size of the 
system obtained from NPT simulations. It is represented in Figure 1, which also 
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Fig. 2. Time-dependence of the collective correlator (q = 0) (solid line) and "self" corre- 
lator (dotted line) for A = A' = (1, 0) and T = 220K, 280K 

contains the experimental result [21]. It is the shrinking of the lattice constant (cf. 
Fig. 1) which leads to an increase of the steric hindrance with decreasing tempera- 
ture and in turn to a slowing down of the orientational dynamics. 

On the other hand, MCT provides an equation of motion for S\y (q, t) which re- 
quires the static correlators S\\> (q, 0) as an input. Because of the crystal's anisotropy, 
these correlators also depend on the direction of q. Therefore the numerical solu- 
tion of the t-dependent MCT equations is rather involved such that we will restrict 
ourselves to the determination of the MCT glass transition temperature T c and the 
normalized nonergodicity parameters 

/xA'(q) ~ lim 4>w{ci,t), (4) 

t^oo 

Md,t) = S X y(^t)/[S XX ( q ,0)Syy(c l ,0)] 1 / 2 

taken at T = T c . They are solutions of an infinite set of nonlinear coupled algebraic 
equations [17]: 

AA'(q)=^AA'(q,{/AA'(q)}) • (5) 
The mode coupling polynomial T\y is related to the memory kernel and depends 
on temperature through the static correlators. The static correlators have been taken 
from the MD-simulation. The numerical solution of Eq. (5) requires a truncation at 
P-max, for which we have chosen £ max = 4. 



3 Results 

As mentioned above the collective correlators (Eq. (3)) are nontrivial for £ > 1 , £' > 
1. Besides these, one can also determine the "self" correlators: 
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Fig. 3. Time-dependence of (a) the collective correlator (q = 0) and (b) "self" correlator 
for A = A' = (1, 0) at T = 300iv" for N = 256 particles (open circles) and N = 2048 
particles (solid line), respectively. 

A N 

sS(*) = ¥ E(W)^'(o)) (6) 

iV 71=1 

iV qel.B.Z. 

8pn,\(t) is the expansion coefficient of 5p n (Q,t) with respect to Y X (Q). Figure 2 
shows the normalized collective correlator at q = and the "self" correlator for 
A = A' = (1, 0) and T = 220K and 280K. First of all we observe a slowing down 
of the a-relaxation by a factor of about twenty by changing the temperature from 
280K to 220K, i.e. by about one fifth. Second, no significant difference occurs 
between the collective and the "self" correlator. The former possesses some oscil- 
lations on the microscopic time scale below 10 ps, in contrast to the self correlator. 
At the lowest temperature 220K both reveal a two-step relaxation process, how- 
ever, with a plateau height very close to one. Since the system size of 256 sites is 
not very huge, we have also performed a simulation at T = 300iv" with 2048 sites. 
The result for A = A' = (1, 0) is presented in Figure 3a (collective) and Figure 3b 
("self") and shows no significant finite size effects. Whether this holds also at the 
lowest temperature T = 220K has not been possible to study because of very long 
equilibration times. 

One of the predictions of MCT is the validity of the t — T superposition principle, 
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Fig. 4. (a) Collective correlator (q = 0) and (b) "self" correlator for A = A' = (1,0) 
versus the rescaled time t/r. r has been determined as the time at which the correlators 
have decayed to l/e. Correlators are shown for T = 220K to 350K in 10K steps. 




t/X 

Fig. 5. Same as Fig.4b, but for A = A' = (3, 0). The inset shows the lowest five tempera- 
tures, only. 

i.e. scaling of t by the a-relaxation time r AA /(q, T) should yield a data collapse 
for all T and t large enough. Taking r AA /(q, T) as the time at which 4>\\>(q_,t) 
has decayed to 1 /e the rescaled data are given in Figure 4a and Figure 4b for the 
collective (q = 0) and "self" correlator, respectively, for A = A' = (1,0). A 
satisfactory collapse is found over a rather huge temperature range. As can be seen 
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Fig. 6. r 1 / 7 versus T for (a) the collective correlator for A = A' = (1, 0) at q = and (b) 
the "self" correlator for A = A' = (1, 0), (2, 0), (2, 2) and (3, 0). 

from Figure 5 this is not true for 0^30 (f). Restriction to 220K <T< 260K leads 
to a satisfactory result for t/r > 1, which, however, excludes the von Schweidler 
regime (see inset of Fig. 5). One of the highly nontrivial predictions of MCT is the 
validity of the power law dependence of r on (T — T c ) : 

t(T) ~ (T - T c y~< , T>T C (7) 

Figure 6b shows r -1 / 7 versus T for the "self" correlators for different A = A' and 
Figure 6a the corresponding result for the collective one at q = and A = A' = 
(1,0). There is a linear behavior for 250i\~ < T < 350i\~. Its extrapolation to zero 
yields 

t md 217 ± 5 K 
For the exponent we find 7aa' = 2.0 ± 0.2 depending on A = A'. 

An interesting behavior has been found for the collective correlators with A = A' 
and i = 2,3. The corresponding results at the Brillouin center are presented in 
Figure 7a and 7b for the lowest temperature T = 220K. Both correlators exhibit a 
high sensitivity on the m-index. For 1 = 2 and I = 3 there is a two-step relaxation 
for m = 0, 2 and m — 0, 1, 3, respectively, with rather high plateau values, partic- 
ularly for 1 = 2. But for t — 2, m — 1 and i = 3, m = 2 there is a rather unusual 
t-dependence showing no typical two-step relaxation with a convex "short" and a 
concave "long" time part joining at an inflection point (the /3-relaxation regime of 
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Fig. 7. Time- and m-dependence of the collective correlator at q = and T = 220 K for 

(a) £ = £' = 2 and (b) £ = £' = 3 

MCT). The "long" time relaxation for t > Aps has a convex curvature on the log- 
arithmic time scale starting at a height of about 0.05 for both cases. If this height 
is interpreted as a plateau height, it has an unusual small value. Surprisingly, it is 
observed that the correlators without two step relaxation decay faster as the tem- 
perature is lowered. 

For (I, m) = (4, 0), (4, 2) and (4, 4) we also have a two step relaxation, while for 
(£, m) = (4, 1) and (4, 3) it is missing, as for (2, 1), (3, 2). The results for i = 4 
yield no new physical insight and are not shown in this work. 

In a final step we have investigated how far these plateau values (which are roughly 
equal to the critical nonergodicity parameters) are reproduced by MCT for plastic 
crystals. First, we have fitted the MD-results by the von Schweidler law in order 
to determine the critical nonergodicity parameters fxx'(<i) from the correlators at 
the lowest simulated temperature T = 220iv", which is close to T C MD . On the other 
hand we have used the static correlators from the simulation to solve Eq. (5) by 
iterations. This allows to locate the MCT-temperature T C MCT , at which the trivial 
solution fxx'(q) = changes discontinuously to fxx'(<l) 7^ 0. 

As a result we have found 

t mct ^ 26 7K 

The corresponding critical nonergodicity parameters from MCT are given together 
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Fig. 8. q-dependence of the critical nonergodicity parameters along the [001], [Oil] and 
[111] directions for A = A' = (1, 0), (1, 1). The result from the MD simulation is given by 
the solid line and from MCT by the dots (the line connecting the dots are a guide for the 

eye) 
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Fig. 9. Same as Fig. 8, but for A = A' = (2, 0), (2, 1), (2, 2) 

with those from the simulation in Figures 8-10. The q-dependence is shown along 
three highly symmetric reciprocal directions. Let us discuss first i — £' — 1 (Fig. 8). 
The numerical result is very close to one for all q. For m = m' = 0, there is a shal- 
low minimum along [001] direction. This minimum also appears in the correspond- 
ing MCT-result. However, / 1 c 010 (q) from MCT is not close to one and deviates 
from the MD-result by about 50 per cent. Furthermore, the former exhibits a much 
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Fig. 10. Same as Fig. 8, but for A = A' = (3, 0), (3, 1), (3, 2), (3, 3) 

stronger q-dependence, whereas the MD-results are practically q-independent with 
exception for m = m' = along the [001] direction. Similar conclusions can be 
drawn for £ = £' = 2 (Fig. 9) and £ = £' = 3 (Fig. 10). There are two features 
which are reproduced by MCT. First, we observe that the MCT-result for / 21 2 i(0) 
(Fig. 9) and /f 2i32 (0) (Fig- 10) is practically zero, in semi-quantitative agreement 
with the corresponding MD-result. This is even true for all q. Second, the hierarchy 
/ 2 c i, 2 i(0) < / 2 C 2 , 22 (0) < /|o, 20 (0) and / 3 c 2i32 (0) < / 3 c li31 (0) < / 3 c 3j33 (0) < / 3 c 0i30 (0) 
is valid for both MCT and MD. 



4 Summary and conclusions 

The main motivation of the present paper has been the test of the validity of the 
MCT predictions for the glassy behavior of molecular crystals. MCT has been re- 
cently extended from liquids to molecular crystals [17]. Since the corresponding 
equations of motion have the same structure as for multicomponent simple liquids, 
the MCT predictions for the liquid systems also hold for molecular crystals. 

In order to check the validity of these predictions we have performed a MD sim- 
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ulation for a simple model of chloroadamantane. In contrast to earlier simulations 
[22,23], we have used a rigid fcc-lattice, however, with an appropriately chosen 
temperature-dependent lattice constant a(T). The decrease of a(T) with decreas- 
ing temperature enhances the steric hindrance, which is responsible for the glassy 
dynamics. Figure 2 clearly demonstrates the validity of this physical picture. 

There are two kinds of tests. On the one hand one can check the validity of the 
two scaling laws predicted by MCT [2,3,4]. Without attempting to calculate the 
corresponding exponents from first principles this represents an important quali- 
tative test. On the other hand a first principle comparison can be made between 
the MD- and MCT-results, as has been done in Refs. [5,6,7,8,9,10,1 1,12,13,14,15]. 
Concerning the test of the scaling laws we have restricted ourselves to the sec- 
ond scaling law which holds for the a-relaxation regime. In this regime, MCT for 
liquids and molecular crystals predicts the t — T superposition principle which 
implies a data-collapse under rescaling time by the a-relaxation time r. Figure 4 
clearly demonstrates the validity of the second scaling law over a large temperature 
range for e.g. the collective and "self" correlators with A = A' = (1, 0). Increas- 
ing £ and £' leads to a shrinking of that temperature range and also of the interval 
of validity for t/r (see Fig. 5). The T-dependence of r is shown in Figure 6. As 
predicted by MCT a power law dependence is found with a glass transition tem- 
perature T C MD = 217K. Close to T C MD deviations from the power law exist, due to 
ergodicity restoring processes. The value T C MD is about 20 per cent below T C MCT , 
consistent with what has been found for liquids [5,6,7,8,9,10,11]. 

The MD-results have shown a very high sensitivity of the relaxational behavior on 
m. For instance the collective correlators for £ = £' = 2 and m = ml = 0, 2 
exhibit a two-step relaxation process with a plateau very close to one, but a very 
peculiar behavior for m — ml — 1 (see Fig.7a). For the latter there is practically 
no plateau. The same holds for £ = £' = 3 and m = ml = 2 (see Fig. 7b). 
A first principle comparison between the MD- and MCT-results (see Figs. 8-10) 
reproduces this unusual m-dependence. We stress that this unusually low plateau 
value is not related to a small corresponding static correlator. However, the critical 
nonergodicity parameters for £ = £ = 2, m = ml = 0, 2 and for I — £' — 3, 
m — ml — 0, 1, 3 from MCT deviate by about 50 per cent from the numerical ones. 
This strong discrepancy can have several reasons. First, the truncation at £ max = 4 
may lead to a significant underestimation of the critical nonergodicity parameters. 
Second, there may be quite a different reason we want to explain now. The potential 
energy of the rigid molecular crystal can be written as follows: 

v(n u . . . , n N ) = v + x: vi(n n ) + \ ]T v 2 (n n , n m ) (8) 

in case of two-body interactions. The isotropic part V is not essential. Besides the 
pair interactions given by V 2 there is a one-particle term Vi, which is a kind of crys- 
tal field, due to the crystal's anisotropy. Vi(Q) will have a finite number of local 
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minima. If the pair interaction V2 would be zero, the molecules at finite temper- 
atures would perform thermally activated jumps between these local minima, as 
can be really seen in Figure 2 of Ref . [22] . Starting in one of the local minima it 
is obvious that the long time limit of, e.g. the "self" correlators, will be positive. 
This happens for all initial conditions with energy below the corresponding energy 
barrier, and has already been investigated for a 4 -model [24]. This type of cage ef- 
fect generated by an anisotropic one-particle potential is not accounted for by MCT 
which describes a two-, three-, etc. particle-cage-effect in a self consistent way. We 
believe that the large numerical values for the critical nonergodicity parameters are 
mainly due to such a one-particle-cage-effect. 

Let us finally comment on the role of phonons, which are absent in the present 
MD-simulation, and also in MCT for molecular crystals [17]. A comparison of, 
e.g. T C MD and 7 from the present simulation with earlier ones for the same model 
for chloroadamantane but on a non-rigid lattice shows that phonons do not have a 
noticeable influence on these quantities and on the long time dynamics. 
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